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Abstract The features of solitary waves observed in 
horizontal monodisperse chain of barely touching beads 
I not only depend on geometrical and material properties 
. of the beads but also on the initial perturbation provided 
' at the edge of the chain. An impact of a large striker on 
, a monodisperse chain, and similarly a sharp decrease of 
' bead radius in a stepped chain, generates a solitary wave 
train containing many single solitary waves ordered by 
I decreasing amplitudes. We find, by simple analytical ar- 
' guments, that the unloading of compression force at the 
, chain edge has a nearly exponential decrease. The char- 
' acteristic time is mainly a function involving the grains' 
, masses and the striker mass. Numerical calculations and 
■ experiments corroborate these findings. 

Keywords one-dimensional granular chain • solitary 
wave • stepped chain • pulses train formation 



1 Introduction 

The problem of mechanical energy transport in an align- 
ment of grains that is held within boundaries has at- 
tracted significant attention in recent years [J; H S S 
0| . It is well established that sound waves may prop- 
agate in such a medium and acoustic speed increases 
with imposed static stress. In the presence of loading, 
linear acoustic waves disperse as a function of time and 
space 0; Q- When the chain is not loaded, i.e., when 
the grains barely touch one another, acoustic speed van- 
ishes: the medium cannot sustain linear acoustic waves. 
In this sonic vacuum limit, nonlinear phenomena pre- 
dominate and any sharp enough perturbation is known 
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to result only in a solitary wave which travels through 
the chain [l|. 

It is well known that mechanical energy transport 
from one grain to another, when the grains are barely 
touching, is a highly nonlinear process. The spherical 
grains interact via the intrinsically nonlinear Hertz po- 
tential, in which the repulsive potential between the grains 
increases with grain-grain overlap via a 5/2 power law. 
This power law is softer than the quadratic (harmonic) 
repulsion for small overlaps but rapidly becomes steeper 
as compression increases. Hence, it is energetically ex- 
pensive for two adjacent grains to sustain contact for 
long. This physical feature is at the heart of the pecu- 
liar nonlinear properties associated with mechanical en- 
ergy transport in granular alignments. It can be shown 
that whenever the repulsive potential has a steeper than 
quadratic power law growth in the overlap (or intergrain 
distance), mechanical energy transfer from one grain to 
the next starts off slowly and then ends abruptly. En- 
ergy is hence transported as a "lump." These alignments 
therefore propagate energy in the form of solitary waves. 
From a geometrical standpoint, the solitary wave in a 
granular alignment of spheres is a highly symmetric ob- 
ject with a well defined center and its velocity depends 
upon its amplitude or on the amount of energy it carries. 
Details of the properties of these solitary waves can be 
found elsewhere [l|; i; 0; H; 1 ■ 

Attention has now turned to mechanical energy prop- 
agation due to initial perturbations that generate signif- 
icant grain compressions at a chain end. When kinetic 
energy is given to a single grain at time t = 0, grain- 
grain interactions force the system to partition this en- 
ergy into part kinetic and part potential energies. The 
energy partitioning between kinetic and potential is con- 
trolled by the nature of the intergrain potential, in this 
case, the Hertz potential. The kinetic energy imparted to 
each edge sphere hence naturally matures into a single 
solitary wave in space and time. When the contact time 
of an impact on the edge grains in a non-loaded align- 
ment becomes sufficiently large, acoustic waves propaga- 
tion facilitates to rapidly distribute the grain compres- 
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sion among multiple contacts. Thus, longer-lived contact 
times at a chain edge can generate comparable grain- 
grain overlaps in several contacts. 

Lazaridi and Nesterenko were the first to demon- 
strate that an initial perturbation provided at the edge of 
a non-loaded granular chain by a large enough striker can 
result in the formation of solitary wave trains, containing 
many single solitary waves ordered by decreasing ampli- 
tudes [1]. These observations have been confirmed and 
extended in a recent simulational work The force 
amplitude of the p*^ single pulse in a solitary wave train 
has been found to decrease approximately exponentially, 
i.e. max ~ exp {—ap). The exponential coefficient 

a was demonstrated to depend only on striker/bead mass 
ratio, even though the attempt to find a simple rela- 
tionship failed In an alternate way of generating 
solitary wave trains, it has also been shown that 
if A'o edge grains are simultaneously set moving from 
their equilibrium positions in the direction of the chain 
at some initial time t = 0, then Nq independent soli- 
tary waves of progressively diminishing amplitudes form 
in the chain. That's exactly what one may expect in a 
Newton's pendulum with several spheres: when two or 
more spheres are used to initiate a s wing , just as many 
spheres swing out at the far end [ill; [12 [H ■ These nu- 
merical findings also agree with observations in a vertical 
column of grains colliding with a rigid wall [l3l | , for which 
the time of collision increases linearly with the number 
of grains contained in the chain due to energy redistri- 
bution of intrinsically dispersive compression waves. In 
another way, several authors have also demonstrated soli- 
tary wave trains may form in gran ular systems contain- 
ing soft and hard materials d; [3; [H; [l3| . 

In this work, we report experimental studies, detailed 
analysis, and numerical simulations, to understand how 
the nature of the initial perturbation provided to the 
edge of a non-loaded granular chain can result in the 
subsequent formation of solitary wave trains. The ex- 
perimental setup under study is a stepped chain, defined 
as two adjacent and non-loaded monodisperse granular 
alignments. A stepped chain is similar, from the point 
of view of the interaction of the last bead of the first 
chain with edge beads of the second smaller chain, to 
a monodisperse chain impacted by a large and massive 
striker [1; |T3; [l^; [13] ■ A stepped chain may also be seen 
as a particular case of a tapered chain (i.e., a chain with 
progressively decreasing bead radius) [l^; [53; [Ml; [H; [1^ . 
Nesterenko et al, then considering a system similar to a 
stepped chain to study the interaction between two sonic 
vacua, reported total energy and momentum transmis- 
sions when a pulse passes through a sharp decrease of 
bead's size, and partial reflection in the reverse case [l3; 
fisj l . Such nontrivial and fascinating behaviors are highly 
related to the way energy and momentum transfer non- 
linearly at the interface between chains. The aim of our 
work is to determine the features of the solitary wave 
train generated when an impulse, initiated at the larger 




Fig. 1 Experimental setup: (top) Schematic view show- 
ing a chain with sensors and data acquisition facilities, (a) 
Monodisperse chain and (b) stepped chain, (c-d) Portion of 
the chain containing a force sensor embedded in a bead (e) 
and ended by rigid and non moving force sensor (f). 

edge beads, crosses the radius mismatch. Starting from 
the conclusions derived in [l^] , we provide here a detailed 
relationship between the shape of the solitary wave train 
and the bead's mass ratio. Making use of the similarity 
mentioned above between the two chain types, our sim- 
ulational studies model a monodisperse chain with large 
striker mass. 

In section ^ we present the experimental appara- 
tus and recall the main elements of a descriptive model 
based on long wavelength approximation of nonlinear 
wave propagation. In section [31 we discuss the exper- 
iments, by means of a qualitative analysis. Finally, in 
section 21 we construct a complementary quantitative 
description and compare it with numerical simulations. 



2 Experimental apparatus and long wavelength 
description of nonlinear wave propagation 

The setup under experimental study is a chain with a 
steep decrease of bead radius, namely a stepped chain. 
All the beads are made of same material (density p. 
Young modulus Y , and Poisson ratio v), but their radii 
Rn and masses m = (4/3)7rpi?^ differ from one chain to 
the other: i?i in the first part and i?2 ^ in the sec- 
ond one. The stepped chain is impacted by a spherical 
striker initially moving at velocity Vs and whose radius 
and mass are Rg and rn^- In our experiment the beads 
and the striker are high carbon chrome hardened ^4/5*/ 
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52100 steel roll bearings with Y = 203 GPa, v = 0.3, and 
p = 7780 kg/m^. Typical elements of our experimental 
setup are shown in Fig. [T] A nonlinear wave is initiated 
in the chain by the striker impact and measurements are 
achieved by using two piezoelectric force sensors. One of 
the sensors is boxed in a rigid and non moving heavy 
piece of metal placed at the end of the chain, and the 
second one is inserted inside a 13 mm radius bead that 
has been cut in two parts. The mass of the bead-sensor 
element has been rectified to match the mass of an orig- 
inal bead 0, and the force measured by the embedded 
sensor is ~ {m-/m)F^ + {m^/m)F-, where -F^ (resp. 
m±) refer to front and back forces (resp. masses) of the 
bead-sensor element Assuming the mass in front of 
the embedded sensor to be negligible (m+ ^ m and ac- 
tually m^/m ~ 0.11) allows to have a satisfactory esti- 
mation of the force exactly at the contact with neighbor 
bead. 

In accordance with pioneering works of Nesterenko [7; 
[3; 0| , the striker impact generates an elastic deformation 
of contact regions between beads, which propagates from 
one to another. The overlap at contact n (i.e. between 
beads n and n—1) is defined as (5„ = {un-i — u„), where 
Un is the position of bead number n. Derived from Hertz 

potential [131 j the force at contact n is F„ = KSf/^ if 
Sn>0 (and if (5„ < 0), where k = R^^^/2^^^e and = 
3(1 — z^^)/4F are constants. More precisely, the dynam- 
ics of the chain is described by applying Newton's second 
law to each of the beads: miln = {Fn — i^„+i). A contin- 
uum limit of this set of equations can be derived under 
the long wavelength approximation, when characteristic 
wavelength A of the deformation is greater than the radii 
of the beads [t^ . Keeping terms of up to the fourth or- 
der of small parameter (2R/X) in Taylor expansions of 
Un±i — u{xn ± 2R,t) versus m„ = u{xn,t) leads to a 
nonlinear equation for the strain ^j}{x,t) = —dxu{x,t), 



^ ~ (2i?)5/2(^j/TO)[tf.3/2 + (2/5)i?2^'l/4(tf'5/4),,],,. (1) 

An exact solution of the last equation is, 



<F^1^m cos-* £. 



vt 
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where one hump of the strain wave 'F, for |^| < tt/2, 
describes a single solitary wave (SW) solution, with a 
compact spatial extent of the order of few times the 
characteristic length A — R^/TO. Thereby, the order of 
magnitude of the neglected terms, (2i?/A)^ ~ 10%, pro- 
vides an estimate of the accuracy of the continuous solu- 
tion, Eq.[21 To complete the long wavelength description, 
one can estimate overlap and force at contact n and ve- 
locity of bead n from the continuous wave description, 
as 5n ^ 2Rx W[x = {2n + l)R,t], F^ = ndl''^ , and 
Vn^vy. if(x = 2Rn,t). 
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Fig. 2 Experiments in stepped chain made with 7 beads of 
radius Ri = 13 mm followed by up to 50 beads of radius 
R2 = 6.5 mm. Force is measured at the end of the chain, and 
the number of beads is progressively increased to follow pulse 
modification as it propagates deeper in the stepped chain, (a) 
Force as a function of time, describing a SW measured at the 
end of a monodisperse chain of 7 beads of radius Ri . This SW 
then generates a SWT when entering the chain of beads with 
radius R2. (b-c) SWT in the second part of the stepped chain, 
respectively at intermediate position (7 large beads followed 
by 25 smaller beads) and at the end of the stepped chain (7 
large beads followed by 50 smaller beads) . (d) Position of local 
maxima of the force as a function of their respective time of 
flights. Wave splitting after passing through the interface is 
clearly shown. 



3 Experimental observations and qualitative 
analysis under weak interaction assumption 

We present experimental observations carried out first 
in a monodisperse chain, and second in a stepped chain. 
The aim of this section is to derive qualitative estimates 
based on the crude assumption that bead interactions are 
described through weakly coupled collisions. We will also 
assume that the interaction between the two parts of the 
stepped chains is governed by the interaction between 
edge beads of each chains. Thus, the last bead of the 
first chain is considered as a striker for the second chain 
with smaller beads. 

As previously demonstrated experimentally and nu- 
merically SHI; Hi], if the initial impact is applied by 
a light striker (m^ <C m), a single solitary wave (SW) 
propagates in the chain. Such a wave is shown in inset 
(a) of Fig. [21 which represents the force measured at the 
end of a short monodisperse chain containing 7 beads of 
radius i?i = 13 mm, and resulting from the impact of 
a striker of radius Rs — 6.5 mm. The profile and the 
velocity both exhibits satisfactory agreements with the- 
oretical counterparts presented in Eq. El [1; S B 01 • 
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A necessary condition for the observation of a single 
SW should be the striker rebound after a single collision: 
the striker won't ever impact the chain again. Applying 
energy and momentum conservation for describing the 
initial striker impact on the first bead of the chain (con- 
sidered uncoupled from the rest of the chain), striker 
velocity K and first bead velocity Vf become after their 
collision, 



Vs and V} 



2m„ 



(3) 



which indicates that striker velocity is negative after im- 
pact if it is lighter than a bead (m^ < to). Then, the first 
bead of the chain accelerates right after striker impact 
and subsequently transmits energy and momentum to 
next nearest neighbor through an almost perfectly nonre- 
flecting collision. From one bead to another, this process 
might generate a wave with compact spatial support in 
the chain, i.e. a single SW. Based on this qualitative de- 
scription, the amplitude of this wave, expressed in terms 
of bead velocity, might thus be roughly. 



2m, 



rus + TO 



Another condition for the observation of a single SW, 
should also be an initial collision duration shorter or 
equal to the duration of the transmitted SW, which is 
approximately two times the characteristic time = 
R^/TO/v. Hence, fraction of initial energy and momen- 
tum will be transmitted through a duration shorter than 
the characteristic time of the chain. The duration of the 
collision between the striker and the first bead (con- 
sidered uncoupled from the rest of the chain) is Tc = 
(2.94/V;^/^)(^/A:)2/5, where n = 1/(to-i + toJ^) is the 
reduced mass and fc-^ = {59/2){R-^ + R-^Y^^ 26]. Us- 
ing first, the relation between v, tPVn and Vm, second the 
relation between Vm and Vs given by Eq. |4l and finally 
introducing the radius ratio rj = R/Rg, the ratio of du- 
rations, 



(5) 



demonstrates that Tc < 2tsw if the striker is lighter than 
a bead (tos < to i.e. ry > 1). 

On the other hand, if the striker is heavier than the 
beads, the impact of the striker leads to the generation 
of a solitary wave train (SWT) Such a SWT, 

resulting in our experiment from the transmission of a 
single SW from the first part of the chain to the second 
part with smaller beads, is shown in insets (6) and (c) 
on Fig. [2] These insets represent the force measured at 
the end of a stepped chain containing 7 beads of radius 
i?i = 13 mm followed respectively by 25 and 50 beads 
of radius i?2 = 6.5 mm. One sees clearly that the in- 
cident SW, shown in inset (a) on Fig. [2l splits into a 
SWT right after passing through the interface between 
the chains. This SWT, considered far from the interface. 
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Fig. 3 Experiments in stepped chain: (a) force versus time 
(incident and interface). One sees first a characteristic com- 
pression of the incident solitary wave followed by longer ex- 
ponential relaxation, (b) Forces at interface normalized by in- 
cident force amplitude versus time normalized by duration of 
incident wave. Matching experimental data to an exponential 
fit for each of the force amplitudes indicates that unloading 



(4) characteristic time is such that r/2r. 
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0.86 ±0.08. 



contains many SWs ordered by decreasing amplitudes 
and behaving individually according to theoretical pre- 
dictions of Eq.O The wave splitting is also clearly under- 
lined in inset id) of Fig. [2l which represents the position 
of few noticeable local maxima in the force signal (i.e. 
the position where force was measured) as a function of 
their respective time of flight. Finally, analyzing the lo- 
cal maxima of the force far from the interface shows that 
the force amplitude of the /c*'* SW in the SWT behaves 
like F„i{k) (x exp(— afc). Matching experimental data 
with an exponential fit, at nine different positions in the 
second part of the chain, one finds ttexp = 0.48 ± 0.03. 

A closer view of what happens at the interface is 
made possible by using the embedded sensor to measure 
the force as a function of time exactly at the interface 
between the chains. This measurement is presented in in- 
set (a) of Fig. [3l which represents the incident force (the 
force measured at the contact before the interface) and 
the force at the interface. One sees that the interface is 
first loaded within a short duration and is then unloaded 
through a relaxation-like process which roughly follows 
an exponential decrease, with characteristic time much 
larger than the loading phase. To obtain a deep insight of 
this process, the force at the interface, measured for dif- 
ferent incident force amplitudes, and normalized by the 
incident force amplitude is plotted in inset (6) of Fig. [3] 
as a function of time normalized by the characteristic 
time of the incident force. This plot shows a collapse of 
the data in a large range of the incident force ampli- 
tude. This demonstrates, first, that the loading duration 
of the interface equals the loading time of the incident 
wave (i.e. half its duration). Second, the maximum force 
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at the interface Fm is proportional to the incident SW's 
force amphtude Fm \ Third, the unloading characteris- 
tic time T of the interface is proportional to the incident 
SW's duration 2ts^(i). In the experimental range of in- 
cident forces, we find F,iiV^™^ = 0.36±0.01 and match- 
ing the force at interface by an exponential fit indicates 
T/2rs^(i) = 0.86 ±0.08. 

As can be seen in Fig. [2] and Fig. [31 the total dura- 
tion of the compression at the interface and of the SWT 
far from the interface is much greater than the duration 
of the incident single SW. This result may be less obvi- 
ous since, from Eq. \5\ the duration of the initial collision 
is of the order of one single transmitted SW duration 
(tc ~ 2tsw if Rs 3> R). In fact, following our qualita- 
tive description, this suggests that few secondary colli- 
sions between a large striker (i.e. the last bead of the 
first chain) and the first bead of the second chain may 
occur. This assumption agrees with Nesterenko's numer- 
ical simulations, demonstrating that a sequence of pulses 
is induced at the interface of two sonic vacua by the de- 
celeration of interfacial particles [l3|. Indeed, one can 
note from Eq. [3] (with rus > m) that the striker is de- 
celerated after the initial collision but continues to move 
forward, while the first bead of the chain is accelerated 
to a higher velocity (V^ > The striker and the first 
bead thus tend to separate due to velocity differences, 
and stress and strain at their contact quickly vanishes. 
Almost entire transfers of momentum and energy from 
the first bead to next nearest neighbor occur during this 
process and the first bead consequently decelerates while 
generating a single SW with an amplitude Vm{^) given 
by Eq. m The amplitude of the first SW, expressed in 
terms of compression force of bead contacts, is such that 
Fm{l) oc Vm^i^)- The first bead thus ends up being 
decelerated, an increase of stress and strain applied by 
the striker consequently occurs, namely a secondary col- 
lision. The velocity of the first bead after this second 
collision is such that Kn(2) = (m^ — m)/{ms + m)Kn(l) 
(see Eq. [21). Following this simple picture, the duration 
between two successive secondary collisions should be of 
the order of the duration of one single transmitted SW. 
The repetition of these secondary collisions (acceleration 
and deceleration of the first bead) leads to the generation 
of successive single SWs whose amplitudes decrease ex- 
ponentially. Force amplitude of the fc*'' transmitted SW 
is thus Fm{k) oc [{iris 
with, 



m)/ [nis + m)]^*^/^ cx exp (— afc) 



a ~ 



if TO, > TO. 



(6) 



Interestingly, this qualitative description of repeated 
secondary collisions predicts the total transmission of 
striker energy and momentum to the second chain, as 
observed by Nesterenko et al 0; [13; [13 j over an infinite 
duration. 

In addition, denoting tk as the instant of the A:*'* sec- 
ondary collision of the striker with the first bead and 
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Fig. 4 Numerical simulations. The figure shows force peaks 
recorded in a bead in a granular chain where the bead is 
sufficiently far from the edge. The force peaks are recorded as 
a function of the peak number. The force decay is found to be 
approximately exponential in this calculation where the study 
has been carried out via direct solution of the Newtonian 
equations for the Hertzian grains. 



introducing t, = {tk — ifc-i) as the elapsed duration be- 
tween two successive secondary collisions of the striker, 
and finally considering that the force at the interface 
decreases approximately exponentially within a charac- 
teristic duration r, one can predict that the unloading 
of the interfacial force versus time is such that. 



Fint{tk) OC exp {-tk/r) with r ~ r*/a. 



(7) 



The elapsed time between two successive secondary 
collisions of the striker can be estimated from the exper- 
imental measurements of t/2t5„(i) and aexp, such that 
t*/2tsj„(i) = a X [t/2t^^(i)] ~ 0.42 ± 0.07. Addition- 
ally, the comparison of the curves' slopes in inset {d) 
of Fig. [21 indicates that Wi(l) ~ V2{k), where Vn{k) is 
the wave velocity of pulse k in chain n. Thus remem- 
bering that Tgyj(^n){k) = RnVW / Vn{k) , the ratio of char- 



acteristic durations r„. 



'(1) 



'(2) 



R1/R2 



2 leads to 



t,/2tsi„(2) — 0.84±0.14. Consequently, the duration be- 
tween two successive secondary collisions of the striker 
equals approximately the duration of one single trans- 
mitted SW, as stated above. However, the elapsed dura- 
tion is slightly amplitude dependent, through v„(fc) (x 

Fm^{k), and slowly increases as the force at interface de- 
creases. Thus, the exponential approximation should be 
considered with care and may remain valid only if con- 
sidered over a limited range of force amplitude at the 
interface, i.e. over a small number of pulses in the SWT. 



4 Dynamical simulations and quantitative 
analysis by means of quasi-particle description 

Dynamical simulations have been performed, in which 
we have allowed strikers of different masses to hit the 
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edge of a sufficiently long chain. The equations of motion 
have been solved via the velocity Verlet algorithm. Dis- 
sipative effects have been ignored in our numerical study 
although such effects can be readily incorporated. One 
finds that sufficiently far from the chain edge the force 
felt by an individual grain can be plotted as a function of 
time for each value of the ratio ms/m. From the conclu- 
sions drawn in last section, it turns out that it appears 
better to plot the force felt as a function of the force 
peak sequence instead of the time. This data is shown in 
Fig.S 



4.1 Single solitary wave energy and momentum: an 
equivalent quasi-particle description 

We present here a complementary model to describe the 
interaction of a striker with an entire chain, based on 
energy and momentum conservation. Before the colli- 
sion with the chain, the striker carries energy and mo- 
mentum, Es = {\/2)msV^ and Qs = rngVs, respec- 
tively. The impact generates a nonlinear wave in the 
chain, and one knows that the hertzian potential en- 
ergy stored at contact n is C/„ ~ (2k/5)i5„ , the kinetic 
energy of a bead is Kn — {l/2)mV^ and its momen- 
tum is Qn = mVn- Assuming that a single solitary wave 
propagates in the chain, one can estimate total energy 
Esw = J2 i^n + Un) and momentum Qstu = I] Qn by 
summing over all the contacts and all the beads. Denot- 
ing Qn/{2R), Kn/{2R) and /7„/(2i?), as momentum, ki- 
netic and potential energies per unit length respectively, 
the total energy and momentum in the chain can also be 
determined by integrating these quantities over the en- 
tire chain length. The kinetic energy Kgw, the potential 
energy Ugw, the total energy Esw, and the momentum 
Qsw of a single SW sufficiently far from the edges are. 
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(8) 



(9) 




{Ws + Ww), (10) 



(11) 



mv 

Qsw = — ^dx^ 

where W„ — J^^^ cos" {£,)d^ is the WaUis integraQ, which 
is such that Wn>2 = [{n — '^)/n]Wn~2 and PVo ~ 7r/2. 

In order to check the accuracy of previous estima- 
tions, one can compare for instance the kinetic to poten- 
tial energy ratio, Ksw/Usw = Ws/Wio = 10/9, to the 
Virial theorem prediction, which states that Ksw/Usw — 

^ http:/ /fr.wikipedia.org/wiki/Integrales_de_Wallis 



5/4 if inter-grains forces derive from the Hertz poten- 
tial. The difference remains weak (less than 10%) and 
is consistent with the long wavelength approximation to 
obtain continuous equation Eq. [T] and solution Eq. O 
This difference could be lowered by calculating higher 
order correction terms of the potential energy from the 
Lagrangian density of Eq. [1] Q . 

Interestingly, the proportionality of kinetic and po- 
tential energies and the spatial compactness of a sin- 
gle SW give an opportunity to describe its dynamics 
as a quasi-particle \VA [Tsl : 0; Defining an equiva- 
lent particle of mass rusw moving at velocity Vsw , such 
that SW momentum and energy are Qsw = Tn^ 
Esw = {l/2)mswVsw, one finds. 



Q 



2 

SW 



f2 X m where f2 



2Esw 

Vsw —Fx Vm where F 



Ws + Wio ' 
Ws + Wio 



(12) 
(13) 



where Vm = v x !fv« is the maximum velocity of beads 
in the chain, and where f2 ~ 1.345 and F ~ 1.385. The 
effective mass agrees with Daraio's and Nesterenko's pre- 
diction [13], who reported a SW equivalent mass about 
1.4 times the mass of a single bead of the chain. The sig- 
nificance of the SW effective mass may be understood by 
considering that the strain wave spatial extent is a few 
beads wide, and msw roughly corresponds to the total 
mass of these beads weighted by the wave envelope (e.g. 
by their own velocities). This effective model may consti- 
tute a simple and helpful tool to predict SW energy and 
momentum transmissions, by using hard sphere binary 
collisions formalism. As a first example, one can predict 
that a striker with mass rus = m and initial velocity Vs 
impacting an uncompressed monodisperse chain made of 
bead with mass m generates a SW with maximum veloc- 
ity Vm/Vs = 2/[{l + Q)F] ~ 0.616. When the SW further 
reaches the end of the chain, the last bead is ejected with 
velocity VjVs = 412/(1 + Qf ~ 0.978. Previous results 
are both in satisfactory agreement with Chatterjee's nu- 
merical results [25], who reported that Vm/Vs — 0.682 
and Vi/Vs — 0.986. As a second example, Nesterenko 
showed that a single SW propagating from a chain of 
light beads to a chain of heavy ones generates sing le re- 
flected and transmitted SWs at the interface fl7|. Ap- 
plying conservation laws on effective quasi-particles' mo- 
mentum and energy gave reflected and transmitted am- 
plitudes estimates in agreement with computer calcula- 
tions [13; 111] - To the best of our knowledge, quantitative 
description of single SW propagating from a chain of 
heavy beads to a chain of light ones has not been yet 
provided. This is presented in following subsection. 



4.2 Solitary wave train energy and momentum 

We estimate here energy and momentum of a SWT. We 
assume first that a SWT is composed of a sufficiently 
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large number of well defined and separated single SWs, 
and second that their amplitudes, expressed in terms of 
the strain, decrease exponentially, ^m{k) oc exp(— 
where ^m(k) denotes the amplitude of the k*^ solitary 
wave, and (3 £ Total energy and momentum of the 
chain are thus, 

oo 

=£;.^„(l)/[l-exp (-5/3/2)], (14) 

k=l 

oc 

Qs^.t = J2 ^^»(^) = Qsw{l)/[l - exp (-5/3/4)]. (15) 

k=l 

Considering that energy and momentum are totally 
transmitted {Eswt = and Qsv^t = Q.s) at long time if 
the mass ratio {nis/m) is large [3;[l3;ll3|5 one finds the 
following relation, 

g2 ^ X 



(16) 



2Es 2E,^t tanh (5/3/8)' 
Using the relation between force and strain, Fm{k) oc 
^^^(fc) cx exp (— afc), so that a = 3/3/2, one finally finds, 
1 + (m/ms)i? 



5„„ 



1 — (m/ms)f^ 



(17) 



where it can be seen that a g implies that pre- 
vious developments and results are physically valid if 
{nis/m) > i7. The comparison between Eq. I17[ simu- 
lations from Fig. [4] and fl0|, and the experiment from 
Fig. [21 is shown in Fig. [5l A satisfactory agreement is 
shown although numerical and experimental SWTs seem 
to contain a finite number of single SWs. A comparison 
between the expression given by Eq. [17] and the approx- 
imate result given by Eq.[6| confirms that the qualitative 
description given in section[3]is meaningful, provided the 
interaction between the striker and the chain is consid- 
ered as an interaction between the striker with a SW 
equivalent quasi-particle of mass rusw = f2 x m, rather 
than with the first bead of mass m. In addition, the result 
presented in Eq. [T7| depicts experimental and numerical 
observations in the sense that the coefficient a does not 
depend on incident SW amplitude ^ld\. 

The analysis above, in addition, sheds some light on 
the shock absorption capabilities of linear chains. Indeed, 
as shown experimentally in Fig. [31 the normalized force 
at the interface of the stepped chain also appeared inde- 
pendent of incident SW amplitude. According to Eqs.fMl- 
[TTI (when > f2 x m), the leading transmitted SW 
(fc = 1) energy and momentum are indeed found propor- 
tional to incident energy and momentum, and respective 
coefficients of proportionality only depend on mass ratio. 



.(1) 



[1 + f2{m/msW 
2fi{m/ms) 



(19) 



1 -|- f2{jn/ms) 

These last results underline the specific and highly 
nonlinear behavior of non-loaded granular chains, since 
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mass ratio m /m 



100 



Fig. 5 Coefficient a, defined such that Fm{k) oc exp (— afc) 
where Fm{k) is the maximum force of the fc*'' transmitted 
SW, when a striker of mass iris collides a monodisper se ch ain 
of beads with mass m. Theory is obtained from Eq. 1171 and 
simulations data comes from Fig. |3] and [T3|. The experimen- 
tal point corresponds to measurements done in the stepped 
chain described in Fig. (2] with ma/m = mi/m2 = 8. Exper- 
imental coefficient is Oexp ~ 0.48 ± 0.03 and theory predicts 
ath = 0.41. 



the larger the striker is, the longer the duration of the 
SWT (i.e. the smaller a is) and the smaller the leading 
transmitted SW's energy and momentum are. This, in 
particular, might be of interest to design efficient blast 
protection systems to prevent from the impact of large 
objects. 



5 Conclusion 

In summary, we have presented a detailed study of the 
formation of solitary wave trains in a non-loaded chain 
of beads. In accordance with existing results, the quali- 
tative description demonstrated that the impact a light 
striker on an alignment of spheres generates a single soli- 
tary wave, while the impact of a large striker induces a 
solitary wave train. The experimental study, the detailed 
theoretical analysis, and the dynamical simulations us- 
ing Newtonian equations revealed that the force peaks 
felt by any grain in the chain decays approximately ex- 
ponentially with the peak sequence, as a solitary wave 
train passed through it. The use of an equivalent quasi- 
particle description, in addition to energy and momen- 
tum conservation, allowed finding a simple relationship 
between the shape of the solitary wave train and the 
striker/bead mass ratio. The prediction is in satisfac- 
tory agreement with simulations and experiment. Pro- 
vided the striker is slightly heavier than a single bead, 
the overall strength of the transmitted wave lowers and 
its total duration rises as the striker mass is increased. 
This nontrivial and unique feature should prove valuable 
when optimizing blast protection capabilities of such a 
system. For instance, given the time and amplitude scale 
of the impulse force, our results help to determine the 
physical dimensions of a bead protection system. 
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